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We present 5j>zteer observations in five wavebands between 3.6 and 24 /im of 
_ an unbiased sample of nine luminous, dusty galaxies selected at 1200 /im by the 

QQ ' MAMBO camera on the IRAM 30-m telescope, a population akin to the well- 

^* • known submillimeter (submm) or 'SCUBA' galaxies (hereafter SMGs). Owing to 

VO ' the coarse resolution of submm/mm instrumentation, SMGs have traditionally 

^ ' been difficult to identify at other wavelengths. We compare our multi-wavelength 

O ■ catalogs to show that the overlap between 24 and 1200 /im must be close to com- 

plete at these flux levels. We find that all (4/4) of the most secure >4a SMGs 
have robust >4cr counterparts at 1.4 GHz, while the fraction drops to 7/9 using 
jh ' all >3cr SMGs. We show that combining mid-infrared (mid-IR) and marginal 

^ ■ (^3o") radio detections provides plausible identifications in the remaining cases. 
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enabling us to identify the complete sample. Accretion onto an obscured central 
engine is betrayed by the shape of the mid-IR continuum emission for several 
sources, confirming Spitzer's potential to weed out active galaxies. We demon- 
strate the power of a 824^x11/ Ss^im versus Sg^ra/ S4,5fj,m color-color plot as a diagnos- 
tic for this purpose. However, we conclude that the majority (~75%) of SMGs 
have rest-frame mid-/far-IR spectral energy distributions (SEDs) commensurate 
with obscured starbursts. Sensitive 24-/im observations are clearly a useful route 
to identify and characterize reliable counterparts to high-redshift far-IR-bright 
galaxies, complementing what is possible via deep radio imaging. 

Subject headings: galaxies: evolution - galaxies: formation 



1. Introduction 

Major progress has been made in our understanding of galaxy formation and evolution 
since the discovery of a significant population of submm/mm-bright, dusty galaxies. The 
first handful of galaxies, discovered behind lensing clusters (Small, Ivison & Blain 1997), 
and in the Hubble Deep Field (Hughes et al. 1998), revealed the diversity of the population. 
SCUBA (Holland et al. 1999), the innovative camera with which SMGs were first seen, 
was commissioned on the James Clerk Maxwell Telescope during 1996-7. Since then, the 
117-bolo meter, 1200-/im MAMBO camera (Bertoldi et al. 2000) has been installed on the 
IRAM 30-m telescope, Spain. Surveys at 1200 /xm are sensitive to the same dusty galaxies 
as SCUBA, yet with a selection function that potentially stretches to higher redshifts and is 
sensitive to lower dust temperatures (Eales et al. 2003; Blain et al. 2004b). 

A moderate fraction of SMGs can be pinpointed by their /iJy-level emission at 1.4 GHz 
(Ivison et al. 2002), allowing us to identify their counterparts in other wavebands. However, 
up to 30% of the population remain undetected in the radio waveband in even the deepest 
maps and the properties of these galaxies remain a mystery. Those SMGs which can be 
localised through their radio emission are usually found in the optical/IR wavebands to be 
faint, morphologically complex systems, often comprising red galaxies with bluer compan- 
ions, as expected for a distant, dust-reddened, interacting star- forming population (Ivison 
et al. 2002; Small et al. 2002; Webb et al. 2003). Recent advances in their study include the 
measurement of the redshift distribution of radio-detected SMGs (Chapman et al. 2003), the 
resolution of most of the submm background by exploiting gravitational lensing (Blain et al. 
1999; Cowie, Barger & Kneib 2002), the identification of X-ray emission from a significant 
fraction of the population (Alexander et al. 2003), the detection — via broad CO lines — of 
coUosal quantities (~10^^ Mq) of molecular gas (Prayer et al. 1998; Neri et al. 2003; Greve et 
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al. 2004a), and the first indications of the strong clustering expected for such massive galax- 
ies (Blain et al. 2004a). These discoveries underline the importance of studying the far-IR 
luminous SMGs for our understanding the formation of massive galaxies at high redshifts. 

For local IR-luminous galaxies, mid-IR data from the Infrared Space Observatory (ISO) 
were used to develop a number of useful diagnostics to differentiate between those powered 
predominantly by starbursts and by active galactic nuclei (AGN) (Genzel & Cesarsky 2000). 
ISO was sensitive to only the most luminous SMGs, so efforts to assess the preponderance of 
AGN in SMGs have so far been reliant on UV/optical spectroscopy, radio observations, and 
the aforementioned X-ray imaging, often with ambiguous results. The fraction of distant, 
far-IR-luminous galaxies with an energetically dominant AGN remains poorly determined, 
although it is expected to be small: even SMGs with unambiguous AGN characteristics 
(Ivison et al. 1998) are thought to have a major starburst contribution to their bolometric 
luminosities (Frayer et al. 1998; Genzel et al. 2003). This finding has been echoed in studies 
of distant quasars and radio galaxies: far-IR-luminous examples are found to be gas-rich and 
their submm emission is often resolved, as expected if their power originates in large part 
from stars (Omont et al. 1996; Papadopoulos et al. 2000; Stevens et al. 2003). 

In this Letter, we use mid-IR continuum imaging of nine MAMBO-selected SMGs in 
the Lockman Hole to investigate the use of mid-IR data to identify the counterparts at other 
wavelengths of this population and demonstrate the power of mid-IR data for determining 
the likely power source (AGN or starburst) in these galaxies. 



2. MAMBO galaxies and the diagnostic potential of Spitzer 

Greve et al. (2004b) reported the results of a deep unbiased MAMBO survey of the 
Lockman Hole, a region mapped earlier by SCUBA for the '8-mJy Survey' (Scott et al. 2002). 
This region was selected for reasons that included its low Galactic cirrus emission and the 
availability of high-quality complementary data at X-ray, optical/IR and radio wavelengths. 
Not surprisingly, the same region was selected for some of the earliest imaging with Spitzer 
(Werner et al. 2004). The data exploited herein were obtained using IRAC (Fazio et al. 
2004) and MIPS (Rieke et al. 2004), reaching noise levels of a = 0.8, 0.8, 2.8, 1.7 and 30/iJy 
at 3.6, 4.5, 5.8, 8.0 and 24 /im with absolute calibration accurate to ±10%. These compare 
with 5-0" confusion limits (20 beams per source) of around 1.0, 1.2 and 56/iJy at 3.6, 8.0 
and 24yum (Vaisanen, Tollestrup & Fazio 2001; Dole et al. 2004). The images cover a 5' x 5' 
region, centered near the SMG, LE 850.01 (Table 1). Data were reduced using standard 
Spitzer data analysis tools, as described by Le Floc'h et al. (2004) and Huang et al. (2004). 
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Our analysis also makes use of deep optical imaging from SuprimeCam on the Subaru 
8-m telescope {R, 3a ~ 27.0, Vega, 2"-radius aperture), XMM-Newton X-ray data (3o" ~ 5 x 
10^^^ ergss"^ cm~^, 0.5-10 keV) and a Very Large Array radio image {3a ~ 15;uJybeam~^, 
1.4 GHz) (Ivisonet al. 2002). 

The Greve et al. (2004b) 1200-;um sample includes nine >3o' MAMBO sources within 
the region covered by the IRAC/MIPS observations. Four are detected at > 4a, all of which 
should be secure. The other five are > 3-a detections, of which up to two may be the result 
of noise fluctuations, Eddington bias (Scott et al. 2002) or confusion. Greve et al. (2004b) 
performed Monte Carlo simulations to determine the reliability, completeness and positional 
accuracy of their source catalog. We list the flux densities, signiflcance and positions for all 
nine sources in Table 1 and note that the positions are accurate to ±4" (95% confidence). 

Fig. 1 shows the SEDs of typical local starburst- and AGN-dominated IR-luminous 
galaxies and several SMGs with well-sampled UV-radio spectra. It is worth reviewing the 
emission mechanisms operating across the observable wavelength range: where our adopted 
sample was selected, at 1200 /im (rest-frame far-IR), we are sensitive to cold dust (~40k) 
created in copious quantities by supernovae (SNe; Dunne et al. (2003)), which re-radiates 
energy absorbed in the UV from hot, young stars. In the radio, at 1.4 GHz, we are again sen- 
sitive to SNe (and hence to recent star formation), via synchrotron radiation from relativistic 
electrons. Contamination via radio-loud AGN is possible and, at /iJy flux levels, is virtually 
impossible to distinguish from pure SNe-related emission. We are also sensitive to AGN 
via X-ray emission from accretion disks and their associated coronae, although determining 
the origin of the X-rays is not always trivial and heavily obscured (A^(H) > lO^^'^cm"^, 
Compton-thick) AGN can evade detection altogether (Alexander et al. 2003). At the short- 
est IR wavelengths accessible to Spitzer (3.6-8.0 /im), we are provided — via photospheric 
emission from stars — with a relatively unobscured measure of stellar luminosity, possibly 
even of stellar mass in more evolved systems out to 2; > 2. At 24 /im (rest-frame ~7/im for 
a galaxy at 2; ~ 2-3) we are sensitive to emission from ~500-K dust in the circumnuclear 
torus of AGN, and to the warmest dust in starbursts. 

Do our Spitzer observations provide sufficient sensitivity to detect dust-enshrouded star- 
bursts and AGN in the distant Universe? In Fig. 1 we illustrate the SEDs of a selection of 
well-studied dusty, luminous galaxies and models (both AGN and starbursts). These demon- 
strate the range of rest-frame mid-IR luminosities for these systems and we select Arp 220 
and Mrk231 to represent extremes of the the starburst- and AGN-dominated sub-classes. 
Fig. 2 shows that even with its steep mid-IR spectrum {a = —2.9, where F^, oc z/"), Arp 220 
would be detected at 3.6 and 24 /xm out to 2; > 1 (see also Blain et al. (2004b)). Mrk231, its 
mid-IR emission boosted an order of magnitude higher than Arp 220 by its AGN, would be 



-5- 



detected out to 2; > 3 at 24 /xm. 

The diagnostic power of the near- and mid-IR bands — the abihty to discriminate 
between starbursts with and without buried, active nuclei — arises from the different physical 
regimes and spectral features probed across the mid-IR waveband (Rigopoulou et al. 2002). 
A steeper slope between rest-frame ~ 3-10 /im is apparent for starbursts than for AGN- 
dominated galaxies. The latter, typified by Mrk231 and SMMJ02399-0136 (Fig. 1), have 
power-law spectra covering rest-frame ~0.2-10/im (Ivison et al. 1998). In contrast, the SED 
of a starburst like Arp220 has a flatter region between 1 and 4/im. Thus for MAMBO 
galaxies at 2; ~ 2-3 (assuming the same median redshift as for SCUBA galaxies — Chapman 
et al. (2003)), the key spectral indices are available in the 3.6-24-/im bands covered by 
IRAC and MIPS. Fig. 3, a 5'24^m/'S'8^m versus 5'8^m/'S'4.5^m color-color plot, shows the tracks 
of Arp220 and Mrk231 as they are moved from the local Universe to z ^ 1, as well as 
the track of a synthetic lO^^-L© starburst from Lagache, Dole & Puget (2003). We can see 
immediately that the starburst and AGN populations are well separated, and that this plot 
has strong diagnostic potential: we can define regions where starburst and active galaxies 
are likely to lie at z = 1-3 and hence use this to classify our SMGs. 



3. Mid-IR counterparts to MAMBO galaxies 

We search for radio and mid-IR counterparts to the MAMBO sources within their nomi- 
nal error circles (illustrated in Fig. 4). We find that all of the four most secure >4o" MAMBO 
sources have robust radio counterparts at 1.4 GHz (as defined by Ivison et al. (2002)). This 
drops slightly to 7/9 robust radio counterparts using all >3a MAMBO sources, somewhat 
higher than the 18/42 detection rate^^ reported by Dannerbauer et al. (2004). The remain- 
ing two have tentative radio identifications at >3cr - but these are not significant enough 
to be reliable. If instead we had used the 24 //m Spitzer imaging to identify counterparts 
we would have identified 8/9 of the >3a MAMBO galaxies with significant 24yum sources. 
Interestingly, all but one of the seven MAMBO sources with robust radio counterparts are 
detected at 24 /xm, the exception being MM J105207. 2+572558. Combining the two identi- 
fication schemes increases our confidence in the identification of the proposed counterparts 
in both wavebands and provides reliable identifications for all nine SMGs in this region. 

Detections at 24 /im (fwhm ~ 6") halve the mm positional uncertainties (to ±2", 95% 



^^Dannerbauer et al. (2004) relied on a 1.4-GHz image which contains several stripe-inducing >60-niJy 
sources, with a noise level 2-3 times higher than the radio data utilized here. Using equivalent data, our 
fraction of radio-detected MAMBO sources would fall to 5/9, consistent with Dannerbauer et al. (2004). 
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confidence) , although the accuracy remains short of that provided by radio imaging (±0.6- 
1.6", 95% confidence — Ivison et al. (2002)). Hence the radio remains the most useful 
waveband for localising the far-IR emission given sufficiently long integrations. 



4. Weeding out AGN 

We next investigate the power of mid-IR observations to identify AGN within the SMG 
population. Only one of our nine MAMBO sources has an X-ray counterpart, MM J105200. 2+572425 
(= LE 850.08), previously identified by Ivison et al. (2002). It could be argued that the com- 
plexity of this immense merger — there are at least seven distinct optical components in 
Fig. 4 — is such that it would be surprising if at least one of the galaxies were not active! 
In the absence of the X-ray detection, could we have identified which components contain 
buried AGN using only the mid-IR data? We could. We list the fiuxes of the three resolved 
mid-IR components in this system separately in Table 1 to check whether they have differ- 
ent SEDs. Fig. 1 shows the individual SEDs, arbitrarily placed at z = 2.5. The component 
near-coincident with the brightest XMM-Newton and radio emission ('b' in Fig. 4; see also 
Ivison et al. (2002)) has a power-law mid-IR spectrum reminiscent of Mrk231. Component 
'a' (previously assumed to be the source of X-rays, based on a less precise centroid from 
ROSAT) was identified as an AGN at z = 0.97 by Lehmann et al. (2001)) and shows a 
weak break in slope at an observed wavelength around 4/im, suggesting a possible mix of 
AGN and star formation power sources. Component 'c', on the other hand, has a strong 
short- wavelength upturn indicative of a starburst. This system allows us to demonstrate 
the diagnostic power of Fig. 3, with component 'b' lying close to the track of Mrk231, well 
separated from the starburst track; component 'a' lies between the AGN and starburst tracks 
and within Az ^ 0.1 of its predicted redshift based on the SED of Mrk231; component 'c' 
is classified as a pure starburst. The main area of possible confusion with this diagram is 
separating high-redshift starbursts from low-redshift AGN; fortunately, as in this case, we 
can expect AGN at 2; ~ 1 to reveal themselves in other wavebands. On the basis of these 
tests, it seems reasonable to expect that mid-IR observations with Spitzer will be able to 
weed out active galaxies within the SMG population. 



The mid-IR data for the other SMGs shown in Figs 1 and 3 display a range of character- 
istics between the extremes of the Arp220 and Mrk231 templates, although the dispersion 
in Fig. 1 must be due partially to the spread in redshift. Over half of the sample display an 
Arp 220-like upturn in the mid-IR — with a fiat SED out to 8 /xm in the observed frame and 
a sharp rise thereafter — showing that active galaxies, although present, do not dominate 
the MAMBO population. This is consistent with conclusions drawn from UV/optical spectra 
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and X-ray imaging (Chapman et al. 2003; Alexander et al. 2003). Using the model/empirical 
tracks to classify the eight MAMBO galaxies with measured 3.6-24-/im SEDs, Fig. 3 would 
indicate six starbursts and two AGN. Thus 75 ± 18% of the SMGs in our small sample have 
mid-IR colors characteristic of high-redshift obscured starbursts. Their mid-IR colors sug- 
gest a median redshift of ~1.4 for our MAMBO sample, somewhat lower than the z = 2.5 
measured for the radio-detected SMGs of Chapman et al. (2003). This may reflect differences 
in the selection functions at 850 and 1200 /im, with the latter selecting a somewhat colder, 
lower redshift and intrinsically lower- luminosity section of the population, or it may simply 
be the result of adopting inappropriate SED templates. 

If we assume the MAMBO sources lie at similar redshifts to the SCUBA population and 
place the MAMBO sources at z ^ 2.5 then 8.0 ;um corresponds to rest-frame K. It is therefore 
instructive to compare the MAMBO galaxies with the i^-corrected spectrum of Arp220 
(Fig. 3). Arp220 has an absolute K' magnitude of —24.4, roughly 1.3 L* (Kim et al. 2002), 
and would have an 8.0-fim. flux of ~2/iJy at z = 2.5. The median 8.0-/im flux of the >4cr 
MAMBO sources is ~20/iJy, implying naively that these galaxies are already ~15L*. The 
more plausible alternatives are that these galaxies lie at lower redshifts than we have assumed 
{z ~ 1-2, as implied by Fig. 3) or that their observed 8.0-fim. luminosity is dominated by 
massive, young stars or hot dust. Determining the exact cause of the apparently immense 
mid-IR luminosities of these SMGs will require spectoscopic observations in the optical and 
mid-IR to confirm their redshifts and the details of their mid-IR spectral properties. 



5. Concluding remarks 

Spitzer observations between 3.6 and 24 /xm of nine galaxies selected at 1200 /im have 
improved our confidence in the identification of secure counterparts, complementing what 
was learned from radio imaging. Robust identifications have been possible for ~90% of the 
MAMBO galaxies using 24-/im imaging, a similar fraction to that achievable from radio 
imaging. Together, radio and mid-IR imaging have yielded plausible counterparts for the 
entire MAMBO sample in this region. We conclude that 3.6-24-/im observations provide a 
useful tool to aid in identifying SMGs. 

Accretion onto an obscured central engine is sometimes betrayed by the shape of the 
mid-IR continuum, confirming Spitzer^s potential to weed out active galaxies. The fraction 
of mm galaxies with energetically-important AGN in our small sample, 25% (or 5-60% at the 
95% confidence level), is similar to that estimated from UV/optical spectroscopy and deep 
X-ray imaging (Chapman et al. 2003; Alexander et al. 2003), but the majority of MAMBO 
galaxies have rest-frame mid-/far-IR SEDs commensurate with obscured starbursts. 



This work is based in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory (JPL), Cahfornia Institute of Technology 
(Caltech) under NASA contract 1407. Support for this work was provided by NASA through 
Contract Number #960785 issued by JPL/Caltech. 
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Table 1. Catalog of MAMBO galaxies in the Spitzer Early Release Observations region 



Source^ 


Si20o(mJy) 


S8S0(mJy)'' 


S24(MJy) 


Ss.oCfJy) 


S5.8{MJy) 


S4.5(MJy) 


Ss.eCpJy) 


Si.4GHz(MJy) 


Comment 


MM J10S204. 1+572658 


3.6±0.6 


9.5±2.8 


249±31 


18.6±1.9 


27.6±3.3 


20.2±2.1 


15.2±1.6 


36±5 


(a); LE8,'')0.14 








— 


18.5±1.9 


22.2±3.0 


14.7±1.5 


9.8±1.0 


72±8 


(b) 


MMJ10S201.S+57244S 


3.4±0.6 


10.5±1.6 


193±31 


10.2±1.7 


11.8±2.8 


8.2±0.9 


3.8±0.8 


73±8 


LB 850.01 


MM J105155. 4+572310 


3.3±0.8 


4.5±1.3 


125±32 


26.4±2.7 


19.9±2.8 


5.8±0.8 


6.1±0.8 


47±5 


LB 850.18 


MM J105200. 2+572425 


2.4±0.6 


5.1±1.3 


534±55 


112±12 


79.0±8.0 


72.6±7.3 


72.0±7.2 


<15 


(a); z = 0.97, LE 850.08 








282±30 


89.5±9.0 


37.2±3.8 


22.3±2.3 


17.1±1.8 


80±8 


(b) 








162±30 


17.8±1.8 


17.2±2.8 


20.7±2.1 


26.7±2.7 


22±6 


(c) 


MM J105157.6 + 572800 


2.2±0.6 


<11 


121±30 


24.5±2.5 


35.7±3.8 


29.8±3.0 


44.5±4.5 


21±6 




MM. J10S207. 2+572558 


1.7±0.5 


<12 


<90 


15. Oil. 7 


14.5±2.8 


17.3±1.8 


12.8±1.3 


23±6 




MM J105203. 6+572612'= 


1.7±0.5 


<10 


166±30 


12.1±1.7 


15.8±2.8 


12. Oil. 2 


7.9±0.8 


19±6 




MM J105216. 0+572506 


1.6±0.5 


6.7±2.1 


578±59 


— 


— 


— 


— 


48±5 


LB 850.29 


MM. J105148. 5+572408 


2.1±0.7 


<7.8 


211±30 


16.7±1.7 


16.8±2.8 


13.4±1.4 


11.9±1.2 


29±6 





^Italics signify sources for which robust counterparts have been identified via radio imaging; bold type indicates a robust identification made possible by the mid-IR, data. 

We quote either a detection or a 3cr limit based on the 850-/J.m noise measured in a 2.6-square-arcmin region centered on the mm source. 
'^The source lies between LE 850.04 and LE850.14; in the SCUBA map it would thus have been affected by the off-beams of both, hence its lack of detection at 850 /J.m. 



-12- 



Rest wavelength /|j,m 
10 




Observed frequency /GHz 

Fig. 1. — Rest-frame SEDs of well-known dusty, luminous galaxies (see Ivison et al. 2000) and 
the 1O^3-L0 synthetic starburst SED from Lagache, Dole & Puget (2003). Mid-IR SEDs for 
the MAMBO galaxies in our sample are shown as dashed lines (vertically offset for clarity), 
normalised to give identical 24-/im luminosity densities and arbitrarily assuming a source 
redshift oi z = 2.5. We distinguish the three components of MM J105200. 2+572421. 9 (a, b, 
c) by plotting them in blue, with a solid line representing the power-law active component, 
'b'. SINGLE COLUMN. 
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Fig. 2. — The variation of flux density with redshift for Arp 220 (pale hnes) and the active 
galaxy, Mrk231 (dark lines). Three bands are shown: 3.6, 24 and 1200 /im, together with 
the 3(7 detection thresholds for our mid-IR/mm observations. As can be seen, our mid-IR 
observations are sufficiently sensitive to identify galaxies with luminosities and SEDs similar 
to either Arp 220 or Mrk231 out to z > 1-3. We assume a cosmology with Qm = 0.3, 
Qa = 0.7 and Ho = TOkms-^ Mpc-^ SINGLE COLUMN. 
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Fig. 3. — A diagnostic color-color diagram showing S'24^m/5'8^m versus 5'8^m/5'4.5^m- Tracks 
for Arp220 (dotted line), a 1O^^-L0 synthetic starburst (solid line), and Mrk231 (an active 
galaxy, dashed line), illustrate how these ratios vary as they are moved from 2; ~ to 2; > 3. 
The diagnostic power of this plot lies in the change in slope of the SED at rest-frame ~3-4 /im 
observed in the spectrum of starbursts. For a submm/mm galaxy at 2; ~ 2.5, this results in 
a pivot around the 8-;um IRAC waveband. Due to their power-law mid-IR continua, AGN 
at z > 1 typically show higher S'8;^m/>S'4.5^m at fixed S24fj,ra/ Ssfj.m than starbursts, providing 
a simple test of the power source. We chose 4.5 /xm for the short-wavelength band because 
it 'comes free' with 8-/im IRAC data; 3.6 /xm would work equally well. Mid-IR data for 
the nine detected MAMBO sources are plotted as points, with the three components of 
MM J105200. 2+572421. 9 plotted separately (identified as a, b, c). These illustrate how 
the active galaxies are separated from the starbursts. In theory, some degree of redshift 
discrimination is also possible via this diagram (see §4). SINGLE COLUMN. 
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Fig. 4. — Multi-wavelength images of 20" x 20" regions centered on the nine MAMBO sources 
within the Spitzer field. For each source we show: left-hand panel, i?-band gray-scale with 
1200-;um contours overlaid at arbitrary levels and, at finer resolution, 1.4-GHz contours at 
2, 3, 4, 5, 6 xa, where a = 5 fiJj beam"^; right-hand panel, 3.6-/ini gray-scale with 24-/im 
contours (at 3, 4, 5, 6, 7 xo"), where the 3.6 and 24-/im images have been aligned. MAMBO 
sources are marked by 8"-diameter circles in the right-hand panels. White squares mark 
tentative identifications based on weak radio emission near-coincident with 24-/im emission. 
FULL-PAGE PLATE. 



